NOTICE 


THIS DOCUMENT HAS BEEN REPRODUCED FROM 
MICROFICHE. ALTHOUGH IT IS RECOGNIZED THAT 
CERTAIN PORTIONS ARE ILLEGIBLE, IT IS BEING RELEASED 
IN THE INTEREST OF MAKING AVAILABLE AS MUCH 
INFORMATION AS POSSIBLE 



NdG- 25269 


Report No. AM-PHY-201-01 

(NASA-CR- 163149) THE ELECTRIC POTENTIAL OF 
PARTICLES IN INTERSTELLAR SPACE RELEASED 
FBOH A NUCLEAR WASTE PAYLOAD Progress 
Report, Jul. 1979 - Feb. 1980 (Alabama ASM 
Univ., Normal.) 45 p HC A03/MF A01 CSCL 03D G3/93 

The Electric Potential Of Particles in Interstellar 


Space Released From A Nuclear Waste Payload 



NASA/George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 35812 


Grant No. NSG-8068 


Unclas 

20929 

l 


February 1980 


Report No, AM-PIIY-201-01 


THE ELECTRIC POTENTIAL OF PARTICLES IN INTERSTELLAR 
SPACE RELEASED FROM A NUCLEAR WASTE PAYLOAD 


by 

Alton C. Williams 
Alabama A. & M. University 
Normal, Alabama 35762 


for 

NASA/George C. Marshall Space Flight Center 
Marshall Space Flight Center, Alabama 35812 


Grant No. NSG-8068 


February 1980 


FOREWORD 


This technical report is the final documentation of work performed 
by Alabama A. & M. University for NASA /Marshall Space Flight Center 
under Grant No. NSG-8068. The period of performance covers July, 1979 
through February, 1980. The results presented in this report are 
intended to assist NASA in its feasibility study of space disposal of 
nuclear waste materials. 
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Report No. AM-PIIY- 201-01 


THE ELECTRIC POTENTIAL OF PARTICLES IN INTERSTELLAR 
SPACE RELEASED FROM A NUCLEAR WASTE PAYLOAD 

1. INTRODUCTION AND SUMMARY 
1.1 Objectives and Scope 

NASA is conducting an initial assessment regarding the desirability 
of disposing of cortuin high-level nuclear wastes in space as an augmen- 
tation of the currently planned national nuclear waste management pro- 
gram. Several space disposal concepts have been considered, including 
solar system escape, injection into the Sun, placement on the Moon, and 
injection into a stable solar orbit [1]. The currently favored concept is 
injection into a circular 0.85 astronomical unit (A.U.) solar orbit, that 
is, about halfway between the Earth and Venus. Orbital calculations 
ndicate that at least for a million years, and probably more, this orbit 
s stable with respect to the Earth and Venus and would not intersect 
the orbits of either one [2], 

The effect of the nuclear waste payload disintegration into small 
particles either by accidental, intentional, or natural means and the 
resulting long-term Earth reentry risk of the particles, however, require 
further analysis. Besides gravitational forces due to the Sun and nearby 
planets, these particles will be subjected to several perturbing forces. 

The most important of these are (1) the Poynting- Robertson drag; (2) 
the bombardment drag arising from the impact of the solar wind; (3) the 
Coulombic drag originating in the momentum transfer from passing solar 
wind components; and (4) the Lorentz force arising from the interaction 
with the interplanetary magnetic field. 

If the grains were electrically neutral, only the first two perturbing 
interactions would have to be taken into account. If they arc charged, 
however, the last two must also bo considered. The relative significance 
of these depends upon the actual electric potential carried by them. It 
is desirable, then, to determine as accurately as possible the potential 
that would exist on such grains. This problem is addressed in this 
report . 


3.2 Summary of Results 

It is found that there are many mechanisms for charging a grain 
in the interplanetary medium. The most important processes are (1) 
capture of solar wind electrons; (2) capture of solar wind protons; (3) 
ejection of electrons through the photoelectric effect due to the solar 




radiation; (4) escape of beta particles from beta emitters in the grain; 
and (5) escape of alpha particles from alpha omittors in the grain. 

If the grains were non-radioactivo (N.R.) , only the first throe 
processes would be present. For such grains, the potentials on them 
would have the following characteristics: Cl) small (£ 12 V); (2) irnle 

pendent of particle size and time; (3) independent of distuneo from Hun 
for distances groutcr than 1 A.U., and (4) for distances less than 1 
A.U., the potentials decrease with distance. Appendix IV contains 
tables for the potentials on grains as a function of the distance from the 
Sun. 


For radioactive grains, the lost two processes would have to bo 
taken into account also. In this cose, the potentials on them huve these 
characteristics: (I) dependent upon particle size, distance from the 

Sun, and time; (2) for spherical grains with radius less than 1 micron 
and located at a distance of 1 A.U. or less from the Sun, the potentials 
are the same as they are ibr the N.R. grains; (3) for grains with radii 
less than or equal to 0.05 micron, the potentials are the same os ibr the 
N.R. grains; (4) for any given distance from the Sun, the potential of 
a groin increases as the radius of the grain increase up to some maximum 
potential (this maximum is equal to the smaller of the potential which 
causes an electrostatic explosion and the equilibrium potential for the 
particular grain); (5) the characteristic time interval over which the 
equilibrium potentials change appreciably is .10 years; and (6) for waste 
which is over 1000 years old, the potentials on grains with a radius less 
than or equal to 10 microns are the same as for N.R. grains as far as 
this solar system is concerned. Numerical results are presented in 
Section 5 where the waste mix is assumed to be PW-4b [3]. 


2. CHARGING RATES 

{ 

2.1 First Order Processes 

In considering the rates at which charge will accumulate on a 
gruin in interplanetary space, various mechanisms of charging must bo 
looked at. The primary or first order processes are defined to bo those 
which can change, the equilibrium potential on a grain by at least 1 V. 
These are (1) the sticking of solar wind electrons; (2) the sticking of 
solar wind ions; (3) the ejection of electrons by the photoelectric effect; 
(4) the escape of alpha particles from alpha emitters; and (5) the escape 
of beta particles from beta emitters, Each of these will now bo looked at 
individually, and the equation which determines the rate of charging for 
each will be given. 


2.1.1 Sticking of Solar Wind Electrons 

The solar wind consists of electrons and ions moving outward from 
the Sun as an aggregate at average speeds which increase as the distance 
from the Sun increases, eventually leveling off at some usymptotie value 
(usually 1 A.U.). The temperature of the solar wind is such that the 
average thermal speed of the electrons is very much greater than the 
mean speed of the wind itself. This means that for the electrons a 
Maxwellian distribution of speeds can be assumed. With this assumption, 
the rate at which units of positive charge will accumulate on a grain of 
radius a is given by [5] 
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e « magnitude of the charge on an electron 
V = potential on the grain. 

The sticking probability, f , for electrons of energies up to about 

V 

15 cV incident upon metallic surfaces lie in the range from 0.5 to 1,0 T5], 

The number density, n o , is about 5 /cm 1 at 1 A.U. and varies inversely 

with both the speed of the solur wind and the square of the solar 
distance [4], Hence, 
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where n Q = 5/cm , = speed of solar wind at r Q = 1 A.U. , and 

V w = speed of solar wind at a distance r. 


2.1.2 Sticking of Solar Wind Ions 


The ions in the solar wind are ubout 90 percent protons. Hence, 
it will be assumed that the solar wind consists of only electrons and 
protons and that it. is electrically neutral overall. The thermal speeds 
of the protons are much smaller than the average speed of the solar 
wind aggregate. It follows , then, that the velocity of the protons is 
the s one as that of the wind. In this case, the rate of gain of positive 
charge on a grain of radius a and potential V is given by [6J 
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where 
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F = fraction of incident protons which remains on the grain and 
P 3 

= number of protons /cm . 


2.1.3 Ejection of Electrons by the Photoelectric Effect 


When radiation from the Sun is incident upon grains, electrons 
are ejected from them through the well known photoelectric effect. 
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The rate at which electrons are ejected depends upon the chemical com- 
position of the grains as well os the potential that they carry. The 
theoretical rate for the rote of gain of positive units of charge for a 
groin of radius a, cutoff wavelength X., and located a distance r from 
the sun is given by [61 1 
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the function which takes into account the fact that the 
ejected electrons will leave the surface with a certain 
energy distribution; 
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V < 0 


2.1.4 Escape of Alpha Particles 


In this section, the rate of gain of positive units of charge due to 
the escape of alpha particles from the alpha emitters in the grain will ho 
calculated. If there are n different such emitters, this rate is given 
by 
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N . " total number of nucleus i present at t - 0 

o,i, a 

X. = disintegration time constant i'or » //-decay of nucleus i 

i $ ti 

F., (a,i) » fraction of the emitted alpha particles from nucleus 1 

1180 which escapes from the grain. 

The function depends upon the potential and radius of 

the groin as well us the range of alpha particles in the grain material. 
This function is determined in Appendix II to lie the following: 
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in the above equations, K, is the average range of the alpha particles 

omitted from nucleus i (see Appendix I) and li. is the kinetic energy 

1 

of those emitted particles, The function 1- is given by 
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Putting equations (2,6) and (2.7) together the rate of increase of 
charge due to alpha emission becomes 
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The number of alpha-emitting nuclei at t ~ 0, N . , can bo written in 

0 ) 1 * 

terms of the mass of nucleus i present at t "0, M, ; the atomic mass of 

1 j 

element i, Aj a ; and Avogadro's number, N^: 
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If the grain hass mass M , and a fraction P. . of this mass is element i, 
then g 1,w 
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If the grain is assumed to be spherical with radius a and has a uniform 
mass density d, then M. can be written, setting M =4/3 ■luftl, as 
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Assuming that this grain was ut one time part of a largor quantity with 

total mass M,„ , with a known distribution of individual nuclei masses, 

m. , then J 
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Equation (2.12) becomes 
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which is the desired equation for the rate of increase of positive charge 
on a grain due to the emission of alpha particles. 


2.1.5 Escape of Beta Particles 

The escape of beta particles from the grain is handled analogously 
to the alpha escape process. However, beta emission is complicated by 
the fact that the kinetic energies of beta particles from any particular 
nucleus form a continuum rather than being discrete as in the alpha 
emitting case. The general equation that determines the rate of increase 
of position units of charge due to beta omission can, nevertheless, be 
written as 
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Whore n 1 is the number of different beta emitters, and the other symbols 
represent either the same or the corresponding quantities for beta 
emitters as they do in the alpha emitting ease. 

The complexity arising from the fact that the kinetic energy spec- 
trum is a continuum is embedded in the function which gives 

the fraction of the beta particles from nucleus i that escapes from the 
groin. This function is determined in Appendix II to be given by the 
following (for a grain of radius a and potential V): 
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whore 


E. = the maximum kinetic energy of the beta spectrum for 
nucleus i: 
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R - R(E) ss the average range of a beta particle with energy E 
in the grain's material (see Appendix 1). 


2.2 Higher Order Processes 

There are various higher order processes which tend to charge 
grains in interstellar space. The most dominant of those result from the 
emission of secondary electrons us churged particles or gamma rays pass 


<) 






through them. The relative significance of the emission of such electrons 
by incident solar wind electrons and protons has already been invosti- 
gated by Wyatt [6], He found that their rate of emission is ubout two 
orders of magnitude less than that of electron ejection by solar quanta. 

Spitzer [71 discussed the emission of secondary electrons by cosmic 
rays and concluded that this rate is negligible compared to the rate of 
sticking of solar wind electrons. Two other mechanisms for secondary 
electron omission must be investigated. Those are (1) the ejection 
caused by emitted alpha and beta particles and (2) the ejection of elec- 
trons by emitted gamma rays. 


2.2.1 Secondary Electron Emission by Alpha and Beta Radiation 

When an alpha particle of speed V collides with an electron at rest, 
it is a straightforward calculation to show that the speed of the electron 
after the collision is given by 
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where 0 is the scattering angle of the electron. Since M/M - 0, this 
becomes e 4 


u ~ 2 V cos 0 


(2.23) 


Hence, the maximum velocity of the electrons is given by 
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This means that the maximum kinetic energy of the scattered electrons is 
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Using E = 6 MeV for the alpha particle energy yields a value of 3.3 keV 

for the maximum electron energy. According to Appendix I, the range 
of these electrons is only 0.027 micron. Hence, even the maximum energy 
electrons will be stopped before escaping from all but the smallest of 
grains. 
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Furthermore, the number of electrons scuttered by on angle 0 is 
proportional to tan 2 0 [8] which means that most of the electrons are 
emitted around 90® with nearly zero energy. In fact, most of the elec- 
trons will be scattered with energies between zero and a few volts, It 
is concluded from this that secondary electron emission by alpha particles 
can be neglected. 

For incident beta particles, the scattered electrons vdll have a 
speed given by U - V cos 0. Because beta particles, in general, are 
much faster than alpha particles, the number of secondary electrons 
emitted by the beta is much smaller. Again, the energy of the majority 
of these electrons is only a lbw volts [13j. Therefore, they will cither 
not traverse to the surface of the grains or not have enough energy to 
escape from the Coulombic attraction of the positively charged grains., 


2.2.2 Secondary Electron Emission by Gamma Radiation 

Gamma particles that are emitted by the various isotopes in the 
grain can interact with the electrons by the photoelectric and the 
Compton scattering processes. These gammas are emitted in conjunction 
with an alpha or a beta particle. Hence, if the number of electrons that 
are ejected by each gamma is not too large, then the isotopes wliich are 
most significant in producing secondary electrons are the same as those 
which are significant for the emission of beta and alpha particles. These 
isotopes are given in the tables in Section 5. 

In determining the relative significance of the secondary electrons, 
it is sufficient to determine the fraction of emitted gamma quanta which 
produces such electrons. This is true because if the gammas originate 
with beta emitters, this fraction is identical with the ratio of the number 
of secondary electrons to the number of primary electrons. Also, if the 
gammas are emitted by alpha emitters this fraction is the ratio of secondary 
electrons to primary alphas which, if small, would mean an even smaller 
secondary to primary electron ratio since the number of alphas emitted 
is significantly smaller than the number of betas. 

For Compton scattering, if N represents the number of electrons 
scattered per unit time, and N q represents the number of incident quanta 

per unit time per unit area then the scattering cross section per electron 
is given by 

a g = N/N Q . (2.26) 

This means that in a length i the number of electrons scattered by the 
Compton effect relative to the number of emitted gammas is given by 


(2.27) 
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where d is the number of electrons per unit volume. The length P, can 
be taken as the radius of the grain. 

An upper limit to the number density of electrons can be obtained 

by assuming that the entire mass of the grain comes from an equal 

number of protons, neutrons, and electrons. If d is the mass density 

of the grain, one obtains for d the following: 
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where M is the mass oi a nucleon (1.67 x 10 gm). Taking d as 
n 3 

6.7 gm/cm and the radius of the grain as 1 micron, the fraction of 
scattered electrons becomes 


f g = (2 x 10 20 ) a g 


(2.29) 


The Compton scattering cross section is a function of the gamma energy 
but is not dependent on the atomic number of the scattering medium [9], 

~25 2 

An upper limit for a„ is about 5 x 10 cm /electron. This gives an 

15 _4 

upper limit for fg as 10 which means that Compton scattered electrons 
are negligible relative to the number of primary electrons. 

For the photoelectric effect, the absorption cross section o per 

atom depends critically upon the energy of the gamma as well as the 
atomic number of the absorbing material [91. In fact a increases as 

the atomic number increases and as the energy decreases. Similar to 
equation (2.27), the number of electrons ejected by the photoelectric 
effect is given by 


fg = yl A l , (2.30) 


where d A is the number of atoms per unit volume. 
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To obtain an upper limit for fg, one can assume that the grain is 

a uniform distribution of a large / element such os lead. Then the 
number density of atoms is given by 



d 



(2.31) 


~22 

where M . is the mass of a lead atom (3.44 x 10 gm). Again, taking 
P D 3 

d as 6.7 gm/cm and l as 1.0 micron, one obtains 


f g = (2 x 10 18 ) o y . (2.32) 

-22 

The range of a in lead for the energies of interest lies between 10 
“20 2 ' -2 

and 10 cm /atom. This gives 10 as an upper limit value for i‘ . 

b 

Hence, as previously discussed since this number is small, the relative 
significance of these secondary electrons is negligible. 


3, EQUILIBRIUM POTENTIALS ON GRAINS 


3.1 General Equation 


The processes which add or subtract positive charge from particles 
in the interplanetary medium have been considered in Section 2. Con- 
centrating on the first order processes, the ones that add positive 
churgo arc the sticking of solar wind protons, the ejection of electrons 
by the photoelectric effect, and the escape of beta particles from the 
beta emitters that may be present. The rate at which each of those 
processes takes place is a function of the potential on the grain. They 
all decrease as the potential becomes more and more positive, and they 
either increase or remain constant as the potential becomes more negative. 


The processes which add negative charge to the grains are the 
sticking of solar wind electrons and the escape of alpha particles from 
the alpha emitters. The rates of these two decrease as the potential 
becomes more negative and either increases or remains constant as the 
potential becomes more positive. 


Consequently, there exists an equilibrium potential at which the 
rate of gain of positive charge and the rate of gain of negative charge 
balance each other. At this potential, the net rate of gain of positive 
units of charge is equal to zero: 



(3.1) 


3.2 Dependence Upon Distance from Sun 

The rates of emission of alpha and beta particles do not depend 
upon the distance from the Sun, but the other processes do. From 
equations (2.1) and (2,2), the rate that electrons will stick to the grain 
is inversely proportional to the square of the solar distance plus any 
dependence upon this distance of the solar wind speed. As was men- 
tioned in Section 2.1 , the solar wind is assumed to be dependent upon 

the solar distance up to some maximum distance r after which it is 

o 

constant. For distances less than 1 A.U., data [4] indicate that the 
average speed of the solar wind is approximately proportional to the 
square root of the distance. Hence, the solar distance (r) dependence 
of the electron capture rate is given by 



(3.2) 
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From equations (2.2) and (2.3), taking the proton number density 
the same as the electron's, it is seen that the rate of capture of solar 
wind protons is inversely proportional to the square of the solar distance. 
Also, as indicated in equation (2.4), the rate of photo-ejection of elec- 
trons is inversely proportional to the square of this distance. 


Hence , 



Substituting equations (3.2), (3.3), and (3.4) into equation (3.1) 

2 

and multiplying by (r/r ) , one obtains the following equation for the 
equilibrium potential : ° 
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(3.5) 


This equation shows that the radioactivity of the grains becomes more 
and more important as the solar distance increases. At very large dis- 
tances, the last three terms become negligible, and the equilibrium 
potential becomes independent of the solar distance. For small distances 
(r < r ), the first two terms which arc duo to the radioactivity of the 

grains become less and. less important. However, even if the radio- 
activity were negligible, the equilibrium potential would still depend upon 

r for r < r Q because of the factor (r /r) n . 
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3.3 Dependence Upon Size of Grain 

Looking at equations (2.1), (2.3), and (2.4), it is seen that the 
electron sticking rate, the proton sticking rate, and the photo-ejection 
rote all ore proportional to the square of the radius a of the grain. 
Thut is, 


where q - e, P, or P.E. Putting this into equation (3.1) and multiplying 
2 

by (a Q /a) > the following equation for the equilibrium potential is obtained: 

(i)' [©„<■>*(")»"] * (S)J, *(s,1 

o o 

+ ( 

The dependence of the radioactive escape rates upon the radius of 
the grain is complex as was shown in sections 2.1.4 and 2.1.5. For 
grain radii that are very much less than the average range of the emitted 
particles, however, these rates are essentially proportional to the cube of 
the radius. In the present case of interest (see Section 5), for a < 1 

g 

micron the radioactive escape rates become proportional to a . This 
means the smaller the grain, the less important the radioactive rates are. 
Conversely, the larger the grains are, the more important the radio- 
activity becomes relative to the other processes. 


dN\ 
dt j 


~ 0 


P.E. a 


(3.7) 


dN_\ 

?“/q 


, v 2 / 

/a \ /' 

vv v 


dN \ 

dt /q 


(3.6) 


3.4 Dependence Upon Time 

For a given grain size at a given distance from the Sun, the 
charging rates due to the sticldng of electrons and protons and to the 
ejection of electrons by the photoelectric effect are independent of time. 
The rates of emission of radioactive particles, however, are dependent 
upon time. This will give the equation for the equilibrium potential a 
non-factorable time dependence and, hence, the equilibrium potential is 
a function of time. 
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According to equations (2.16) und (2.17), the time dependence of 
the rates of charging due to the emission of beta or alpha particles is 
given by 



This shows that each emitter i contributes to the charging rate with a 
time dependent term D.(t) which decreases exponentially with a time 
constant . 1 


1 half-life 
Aj ~ Ln (2) 


(3.9) 


Hence, at times near t - 0, the shorter half-life nuclei could contribute 
appreciably, but as time goes on. only nuclei with relative long half lives 
would contribute. 


Equation (3.8), however, only takes into account the fact that 
parent isotopes are being depleted because they are decaying away. It 
does not include the fact that daughter nuclei are building up due to 
this decay of the parent. If the daughter nucleus is radioactive, it must 
also be included in the terms in equation (3.8). The significance of this 
can be really pronounced if the parent is short-lived and the daughter 
is long-lived. To take this into account quantitatively , one must look at 
the differential equation which governs the number of daughter nuclei 
present at an arbitrary time. 


Let Np and be the number of parent nuclei and daughter nuclei 

present at time t, respectively. Also, let N and N dQ be these same 

numbers at the initial time t = 0. If the parent decays with a disintegra- 
tion constant A^ and the daughter decays with the constant A ( j, then the 

equations that give N p and are the following: 


N = N e p (3.10) 

P po 


d N { 
J dt 


■A, N. + A N 
d d p p 


(3.11) 
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1 


Solving equation (3.il) for N . with the initial condition thnt N, = N, 
ut t * 0 yields a a ao 


A 

t “* X ,, t\ 


N ( G ■ G 1 

a A ,i - A n P° \ / 


* N do c 


V 


(3.12) 


The rate of decay of the daughter nucleus is tho first term In equation 
(3.11), or 


(L**\ 
\ dt / 


decay 




(3.13) 


Using equation (3.12) for N. in equation (3.13) gives the following 
equation : 


( 5 ) 


decay 


“ “1 , N i e 
d do 


V 


1 + / ._J2. J\ U>°. 

1 X , “ A N, 
d p/ do 


( A. \ -At 

\ A N e p 

\l “ k 9 ) p P° 

(3.14) 


This equation is of the form 




= A 


decay 


( 't^)° + B ( d A>)° 

V dt Loay \ dt /d 


(3.15) 


decay 


where (dN/dt) 


is the decay rutc when buildup is not taken into 


decay 

account. Hence, the net effect of the buildup of the daughter nucleus 
is thut the constant in equation (3,8) for this nucleus gets multiplied 

bv the factor A, and *he constant D. for the parent nucleus is multiplied 
by B . That is , 1 


Dj (parent) •+ B D,. (parent) 


(daughter) A D,, (daughter) 


( 3 . 16 ) 

( 3 . 17 ) 
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whore 


A 


end 




N. 

N 


M 

do 



(3.18) 


(3.19) 


Note that the substitutions (3.16) and (3.17) could possibly include 
the following cases: (1) the daughter nucleus may not have boon present 

before the buildup began, and (2) the parent nucleus may be a different 
type of emitter than the daughter, but it will still be included in equu u 
tion (3.8) as if it were the same. 

According to equation (2.13) the number of nuclei present at t ■■ 0 
can bo written in terms of the initlul mass M , the atomic mass A, and 
Avogadro's number as 



(3.20) 


Hence, the constant A cun be written os 



M „o A d 

M do V 


(3.21) 


In summary, the equilibrium potential is time dependent because of 
the rudioaetivity of the groins. The rotes of charging due to this radio- 
activity is given by equation (3.8). The constant C. in this equation is 

obtained for alpha and beta emitters by computing it with equations (2.16) 
and (2.17), respectively. When buildup of any particular daughter 
nucleus is taken into account, the constants C. for the daughter and its 

parent must be multiplied by the factors A and B, respectively. Some- 
times, this may require that either the parent or the daughter be added 
to the equation. In the case of the parent , it may have to be added 
even if it is not the same type of emitter ns the daughter. Jt should 
also be stated that, in general, chains longer than two nuclei would have 
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to bo considered, Those longer chains, however, were ruled out for the 
specific application of this work by checking the effects of other mem hors 
of the chain numerically. 


3.5 Maximum Potential Without Electrostatic Explosion 

For a given grain size and composition, thoro exists a maximum 
electrical potential which it can withstand before it breaks up. The mag- 
nitude of this potential depends upon the tensile strength of the grain. 

To calculate the maximum potential, it will bo assumed that ut the 
critical point the total charge Q is divided into two equal parts of Q/2 
each. Also the two ports are assumed to be separated by one radius 
length a. The Conlombic force of repulsion between the two parts is 
then given by (CGS units) 


F ss - Q-a- . (3.22) 

c 4o 


The force that is tending to hold tho two parts together is the tensile 
strength T times the cross sectional area: 


P t = T irn 


(3.23) 


At the point of breakup, these two forces arc equal: 


P s= TV 
C "t 


(3.24) 


or 


4a 2 


= T ?ra 


(3.25) 


This says that the maximum potential then is given by 


V = S. = 2a (it T)* 

max a 


(3.2G) 
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8 2 

The tensile strengths for most metals He in tho range 10 dyno/cm 

10 2 9 2 

to 10 dyne/cm [10], lienee, if a typical value of 2 * 10 dyno/eo) 

is taken, then 

V max 3 ( 1,585 x ,q5) w * (3.27) 

If a is expressed in microns and V mQX is in V, then 


V max e (<1655)a v • (3,28) 

* 

Since the maximum potential depends upon the square root of the tensile 
strength, even if it changed by a factor of 4 or 5, tho maximum given in 
equation (3,28) would still bo larger than tho potential expected on groins 
(see Section 5). 
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4, RELATIVE SIGNIFICANCE OF LORENT55 AND 
GRAVITATIONAL FORCES 


It is of interest to determine) the relative significance of tlie Lorentz 
and gravitational forces on a grain in interplanetary space ns a function 
of the grain's size and voltage. For a grain of radius a, density d, and 
charge Q , the magnitudes of the solar gravitational and Loroniz forces 
are given by (COS units) 



(4.1) 


(4.2) 


respectively. Here, r is the distance from the Sun, M g is the Sun's 

mass, v is the solar wind's speed, B is the component of the inter- 
planetary magnetic field which is perpendicular to the direction of the 

IQ 

solar wind velocity (radial), and c is the speed of light (3 x 10 cm/soc). 

Both B and v depend upon the solar distance r, but B can bo 
written in terms of v as [11] 


B 

B a -2- , (4.3) 


where B„ is a constant. There are field fluctuations as u function of 
o 

time, but for comparison purposes typical values can be used. Hence, 
at the orbit of the Earth, a typical value for the magnetic field is 3 * 

-5 

10 gauss and the speed of the solar wind is 450 km /sec. This gives 

I 

gauss . (4.4) 

and 


B 




3 x 10 


•5 V o 


>) 


i 
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c 


(4.5) 


F. 


. 3 x 10 


Q v 


o 


«P 
r ’ 


i 


Putting in v Q « 450 km/soc * 450 * 10 5 cm /sec, this becomes 



1350 


Q l ‘o 


(4.0) 


whore r = 1 A ,IJ . 

Taking the ratio of equation (4,0) to equation (4.1), one obtains 


F 

ST “ Cl. 81 v IQ" 8 ) S. , (4.7) 

l G n a cl 


O 

where r is measured in A.U., and u und d are in cm and gm/cm u , 
respectively. Since Q/n is the voltage V on the grain, equation (4.7) 
con be written as 



(0.006) V r 



(4.8) 


where V is measured in V; r in A.U.; a in microns; and d in gm/cm . 

g 

Taking the reference value of 6.7 gm/cm us d, this ratio becomes [3] 


I? 

F, 


(8.96 x 10" 4 ) 

a 1 


(4.9) 


The Loren tz force can be said to be significant when Fj /l’ ( , is 

greater than or equal to 0.1. Using this as a criterion, one can line! the 
maximum potential that a grain of d given radius and solar distance can 
have without perturbing the gravitational force significantly. This maxi- 
mum is given for various radii and solar distances in Table 1. The 
equation used, in general, is 


V 


max 


111,61 ~ 


(4.10) 
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TAHLli 1. MAXIMUM POTENTIAL THAT GRAINS CAN HAVE WITHOUT 
PERTURHING THE SOLAR GRAVITATIONAL FORCE SIGNIFICANTLY 
(r = solar distance; u = radius of grain) 


r 


(A.U.) 

a (microns) 

V (V) 

max ' 

0.1 

0.01 

0.11 


0.10 

11.16 


1.00 

1116.10 

0.5 

0.01 

0.02 


0.10 

2.23 


1.00 

223.22 

1.0 

0.01 

0.01 


0.10 

1.12 


1.00 

111.61 

5.0 

0.01 

0.002 


0.10 

0.22 


1.00 

22.32 


10.00 

2232.20 

10.0 

0.01 

0.001 


0.10 

0.11 


1.00 

11.16 


10.00 

1116.10 


5. NUMERICAL RESULTS WITH PW-4b REFERENCE WASTE 


5.1 Significant Beta and Alpha Emitters 

As was discussed in Section 3, the equilibrium potentials on the 
grains is found by solving equation (3,1). The radioactive decay rates 
in this equation, however, are time dependent as described in equation 
(3.8). In the PW-4b waste mix, there are many radioactive isotopes, 
but only a lew of them contribute significantly to the decay rates at any 
given time. 

Neglecting the buildup of daughter nuclei, for the moment, the 
isotopes which contribute significantly are shown in Tables 2 through 7. 

TABLE 2. BETA EMITTERS FOR 0 < t < 200 YEARS 


Isotope 

E (Mi 

max 

c 90 
Sr 

0.546 

Y 

2.270 

„ 134 
Cs 

0.658 

,, 137 
Cs 

0.514 

r, 147 
Pm 

0.225 

„ 154 
Eu 

0.580 


TABLE 3. BETA EMITTERS FOR 200 < t < 500 YEARS 


TABLE 4. 


Isotope 

E (MeV) 

max 

Sr 90 

0.546 

Cs 137 

0.514 

C 151 

Sm 

0.076 

ETA EMITTERS 

FOR 500 < t 

Isotope 

E (MeV) 

max 

Q 90 
Sr 

0.546 

Cs 137 

0.514 

s ™ 151 

0.076 

ri 93 
Zr 

0.063 

m 99 
Tc 

0.292 

n l 35 
Cs 

0.210 
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TABLE 5. BETA EMITTERS FOR 700 < t ' 7 x 10 6 YEARS 


Isotope 

E (M 

max 

Sm 151 

0.076 

Zr 93 

0.063 

Tc" 

0.292 

Cs 135 

0.210 

Pd 107 

0.035 


TABLE 6. ALPHA EMITTERS FOR 0 < t < 1000 YEARS 


Isotope 


Am 

Cm 


241 

244 


E (MeV) 

5.49 

5.81 


TABLE 7. ALPHA EMITTERS FOR t > 1000 YEARS 


Isotope 

E (MeV) 

Am 241 

5.49 

Np 237 

4.52 

Pu 239 

5.16 

„ 240 
Pu 

5.17 

a 243 
Am 

5.28 


Taking 1 into account the buildup of daughter nuclei, it is found 
that the decay schemes for which either the factor A is significantly 
different from unity or the factor B is significantly different from zero 
[see equation (3.15)] are the following: 


Sr 90 8~ Y 90 3 + 

(5.1) 

Am 2 ^ -> a Np 337 ->• a 

(5.2) 

Np 237 -> a + Pa 233 + R - 

(5.3) 
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For the reaction in equation (5.1), A = -0.25 and B = 1.0. This 

00 

means that for the times when Y is important, its contribution is 

90 

multiplied by -0.25, and when Sr is important, another term will be 

90 

added to the beta emitters which has the same properties as Sr , but 

90 

the energy of the emitted particles is that of Y . 

For the reaction in equation (5.2), A = 1.167 and B ~ 0. Hence, 
237 

at times when Np is important, its scaling* factor D.(t) in equation 
(3.8) should be multiplied by 1.167. 

For the reaction in equation (5.3), A = -0.86 and B =1.0. There- 
233 

fore, when Pa is important, its scaling factor is to be multiplied by 
237 

-0.86. Also, Np must be added to the beta emitters although it is an 

an alpha emitter. When it is added, it will retain the same scaling factor, 
but the energy spectrum of the emitted particles is to be taken as the 

233 

same as that of Pa . It does not become significant, however, until 
after 500 years. 


5.2 Equilibrium Potentials on the Grains 

Solving equation (3.1) numerically, the equilibrium potentials on 
grains as a function of time, grain size, and solar distance were deter- 
mined. These results are tabulated in Appendix IV. By comparing the 
potentials found there with the maximum potential in Table 1 or equation 
(4.10), one can determine for which size grains one has to be concerned 
about the charge on them for a given solar distance and time. For 
example, assume one wants to determine lor what times and distances 
that the charge on a 5 micron radius grain become important. Using- 
equation (4.10), it is found that the maximum potential that the grain 
can have without becoming significant is 558.05 V at r = 5 A.U. and 
279.03 V at r = 10 A.U. From Appendix IV, it is seen that at r = 5 
A.U., the potential on this size grain is less than the needed maximum 
for all times, and at 10 A.U. the potential is significant for all times less 
than 20 years. 

It can also be seen in this way that for solar distances less than 
or equal to 2 A.U. and for all times, the maximum size grain with which 
one has to be concerned has a radius of 0.5 micron. 
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APPENDIX I 


RANGE OF ALPHA AND BETA PARTICLES 


A. Beta Particles 

The range-energy relation for electrons is not strongly dependent 
upon the atomic number of the stopping material [12]. This means that 
the ranges for all metallic type materials are essentially the same, and 
hence one can use the practical range of aluminum as found experi- 
mentally [13]; 



0.412 

d 


,n 


cm 


E < 2.5 MeV 


(.530 E 


. 106) cm 
d 


2.5 MeV < E < 20 MeV 


( 1 . 1 ) 


where, E is the kinetic energy of the beta particle, d is the density in 
3 

gm/cm of the stopping material, and n = 1.265 — 0.0954 Ln E. 


B . Alpha Particles 

The ranges of heavy particles such as alpha particles do depend 
upon the atomic number of the stopping material.. The stopping material 
that is of interest in this work is heterogeneous; hence, an average value 
for the atomic number must be obtained. The particular manner in which 
this averaging is carried out is really not very important since it is found 
that the contribution due to the escape of alpha particles is actually a 
second order effect (see Section 5). 

If one simply calculates the average atomic number of the PW-4b 
mix [3] as it exists in oxide form, a value around that of lead is obtained. 
Hence, for this work the ranges of the alpha particles will be those for 
lead. 

In general, the range of alpha particles in a material is related to 
the range of protons in that material by [12] 

v e) = E p(!) • , (i - 2 > 


The energies E of the alpha particles that are of interest all lie in the 
range from 4 to 8 MeV. Hence, the corresponding proton energies are 
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in the range from 1 to 2 MeV. The range of protons in lead at the 

2 2 

energies 1 and 2 MeV are 7,90 mg/cm and 2.505 mg/cm , respectively 
[12], It is reasonable to assume that the average range of protons 
between 1 and 2 MeV as a function of the energy is a straight line. 

This gives 

It (E') B (17.15E' - 9.25) mg/cm 2 ; 1 MeV < E' < 2 MoV. (1.3) 
Equation (1.2) then becomes 


R (E) = (4.288 E - 9.25) mg/cm 2 ; 4 MeV < E < 8 MeV . 

Ut 


In terms of cm, the range of alpha particles can be wiatten using the 
density d of the material as 


R (E) s 

a 


(4.288 ^ 10“ 3 E - 9.25 x 10~ 3 ) cm 


APPENDIX II 


FRACTION OF ALPHA AND BETA PARTICLES TO 
ESCAPE FROM A FINITE SIZE GRAIN 


Consider the situation where alpha or beta particles are emitted 
uniformly throughout a spherical volume of radius a. Each particle is 
assumed to have initial energy E q and average range R - R(E q ) • 'I'bon 

if a particle travels a distance X (< 2 a) to reach the surface, its energy 
at the surface is given by 


E = E q (1 - X/R) 


(II. 1) 


The fraction of pai’ticles f(X) dX which travel a distance in the 
range (X, X+dX) to reach the surface is given by the fraction of the 
total volume that contains points with distances to the surface in this 
range. It is a straightforward calculation to show that this volume is 
given by 


V(X, 

X+dX) 

dX 

= 

2tt X 2 (1 - X /2a) dX 

(11.2) 

and, hence 

f(X) 

dX = 

3 

2 

X 2 

T 

a 

(1 - X /2a) dX . 

(II. 3) 


The fraction of particles f(E) dE with energy in the range 
(E, E+dE) is the same as the fraction that travels a distance in the 
range (X, X + dX). That is, 


f(E) dE = f(X) dX 


(II. 4) 


Using equation (II. 1) to obtain 


X = R( 1 - E/E ) 


(II. 5) 


and 


(II. 6) 


dX 

dF 


-R/E q , 


equation (II. 4) gives 


f'(E) <1K = 12 


for ull E such that 


(E„ - E)‘ 


o 


E 


(ft)' (' 


i - a - e/e q ) 2 


h) aB 


(II. 7) 


(1 - E/E q ) R /(2a) < 1 


or 


E >12 (1-2 a/R) 

— o 


(II. 8) 


Condition (II. 8) takes into account the fact that some particles may lose 
all their energy before they reach the surface. 

If there is some minimum energy E that the particles need to 

have at the surface to escape, then the fraction of the emitted particles 
with initial energy E q to escape is given by 


F (E 


min' 


V 



E . 
min 


f'(E) 


dE 


(II. 9) 


Using equation (II. 7) for f(E), this fraction is easily found to be 


F (E . , E ) 
v min o 


4 O 


e . n< 

mm 


0 

1 


V 3 (if [ 


1 - - (1 - E - 
4 min 


E . > E 

mm - o 


in'V (I)] 


E . < E (1 - 2n/R) 

min - o 

(II. 10) 
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The fuct that l ? < E min »E ) 85 1 when E jn < K (1 2u/R) comes from 

condition (II. 8) which indicates that nil the particles reach the surface 
with a certain minimum energy. If that minimum is negative, then all 
the particles do not reach the surface. 

Since the initial energies of beta particles omitted from a nucleus 
form a continuum with number density N(E q ) (see Appendix III), whereas 

the emitted alpha particles are essentially monoencrgetic with an energy 
1S Q , then it is necessary to look at alpha particles and beta particles 

separately. Also whether particles that reach the surface will escape or 
not depends upon the sign and the magnitude of the electric potential V 
on the surface. Hence, it is necessary to examine four different enses: 

1) Alpha Particles; V > 0. 

In this cose, the minimum energy needed for an alpha particle 
is zero. Hence the fraction to escape is given by 


F Esc ^ F(0 ’ V ‘ (II,11) 


2) Alpha Particles; V < 0. 

When the potential is negative, the alpha particles that reach 
the surface need a minimum energy of 2 e |V | to overcome the Coulombic 
attraction. The fraction to escape is then 


F Esc = F ( 2e M>V • (II,12) 


3) Beta Particles; V < 0. 


In general, beta particles are emitted with an energy distribu- 
tion N(E q ) where the energies range from zero up to a maximum E max * 

In general, then, the fraction of beta particles to escape is given by 


E 


r max 

J F(E . , E ) N(E ) d E 
J min o o o 


p f P R ) 

A Esc' min’ max^ 


o 


J 


E 


. (11.13) 


max 


N(E o ) d E 0 
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Bor V <_ 0, the minimum energy needed is zero and, therefore, the 
fraction of beta particles to escape is given by 


"lisc <'• ‘W> 


4) Beta Particlos; V * 0. 


In this case, the beta purticlos need a minimum energy of eV 
to overcome the Coulombie attraction. Hence, the frnetion to escape in 
this case is 



! 

f 



APPENDIX III 

ENERGY SPECTRUM OF JJKTA PARTICLES 


A. Momentum Spectrum 

The momentum spectrum of beta particles emitted by a nucleus with 
atomic number 55+1 and mess number A is given by [14] 


N’(P) dP “ C(Z,A) F(Z,A,E) P 2 (E MAX - E) 2 dP , (III. I) 
where: 

N'(P) dP k number of particles emitted with momentum in the 
range (P, P+dP) 

C(Z,A) « a constant for each nucleus (depends upon the length 
of time the nucleus has been emitting particles 

E MAX ” maximum kinetic energy of the beta particles 

F(Z,A,E) « the Coulomb factor which to 1 percent accuracy is 
given by 


F(Z,A,E) 


(2 P f / mo )* B ~ 2 (S 2 + „ 2 ) 8 -* 

( ( 2S ) ! ) ^ 


x Exp (2i|)T1 - 2S + S/(G(S 2 + n 2 ))) 


where: 


S 

= (1 - (Zo 2 /MC) 2 ) 4 

V 

= R/HC; R * nuclear radius 

n 

25 e 2 

= ; v = beta particle's speed 

«l> 

= Tan” 1 (S/n) . 


(III. 2) 

(III. 3) 

(111. 4) 

(111. 5) 

(111. 6) 
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K. Transformation to an Energy Spectrum 

To transfer equation (III.l) to an energy spectrum, one uses the 
fact that the number of particles with energy in the range (E, E+dE) is 
equal to the number of particles with momentum in the range CP, P+dP). 
That is, 


N (E) dE « N' (P) dP . (HI. 7) 

Hence , 

N (E) = N' (P) * N' (P) (dli/dP) 1 . (III. 8) 


The relativistic kincmuticnl equations which are of interest here ore the 
following: 


E “ me 2 (7-1) 

(III* 9) 

P - m yv 

(III. 10) 

/i 2 . 2v-i 

Y - (1 - v /c ) 

(III. 11) 

These can be used to give 

N (E) - C (Z.A) m 2 e F(Z.A.E) (1 + E/mc 2 ) 

[(1 + E/mc 2 ) 2 

- lli(1 Wx- E)2 • 

(III. 12) 

If one defines 

B » (1 + E/mc 2 ) 

(III. 13) 

and 

2 

a - e /be 

(III. 14) 
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then the quantities which define the Coulomb factor I*’ can bo written os 



. J n « 

djT'jJT 

(III. 15) 


s a - z 2 A 1 

(III. 10) 


<p B Tan** 1 (S/n) . 

an. i7> 

Also , 

if a new overall constant is defined by 



C’(Z,A) * - (2 i p) e m c C (Z,A) 

( (2S) ! ) ^ 

(III. 18) 

then 

the energy spectrum is given by 


N(E) 

dE = C' (Z,A) B(B 2 - 1) S “* (S 2 + n 2 ) S ~* t(B max ~ 

E 0 )/mc 2 ) 


X Exp (2 <f>n + S/CCS 2 + H 2 )] dK . 

(HI. 10) 


This is the desired equation for the energy spectrum. 


APPENDIX IV 


EQUILIBRIUM POTENTIAL ON CRAINS 


I. Notation 

a = radius of grain (microns) 
v = distance from the Sun (A. II.) 
t s= ago of reprocessed waste (yr) 

V e(J s Equilibrium Potential (V) 

II. Reference Values Used 

1) Waste form os given try Battolle Laboratories Progress Report 

(September 14, 1979) pg. 3. 

2) Waste Mix — Commercial PW-4b. 

3) Waste reprocessed after 10 years from reactor, 

III. Equilibrium Potential 


A) Submieron Particles (a 1,0) 

1) r * 1 . 0 

Here, the radioactivity of the particles is not significant 
enough to affect the equilibrium potentials. This means that the voltages 
are time independent and ore the same as they are for non-radioactive 
particles : 


r 

..m 

1.00 

12.0 

0.75 

10.5 

0.50 

9.5 

0.25 

8.0 

0.10 

7.0 


2) r > 1.0 

The relative importance of the radioactivity of the particles 
becomes more significant as r increases (also as a increases). 
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? 

fl 

V 

'( 





(i) a = 1.0 


r 

t 

V 

...<a 

2.0 

0 

14 


to 

14 


100 

12 


N.R, 

12 

5.0 

0 

32 


10 

35 


100 

12 


N.R. 

12 

10.0 

0 

175 


10 

95 


100 

17 


1000 

12 


N.R. 

12 

50.0 

0 

5,000 


10 

2,700 


100 

300 


1000 

12 


N.R. 

12 


(ii) a = 0.1 


r 


t 


V 

_eg. 


2.0 0 12 

N.R. 12 

5.0 0 13 

10 12 

N.R. 12 

10.0 0 23 

10 16 

100 12 

N.R. 12 

50.0 0 450 

10 250 

100 25 

1000 12 

N.R. 12 




(iii) u < 0,05 

For particles of this size , the radioactivity is not 
significant for uny reasonable distance from the Sun. (r ^ 100) 

1J) Particles larger than 1 micron (a > 1.0) 

1) r < 1.0, a = 10.0 


r 

t 

js a 

1.00 

0 

19.0 


10 

15.5 


100 

12.0 


N.R. 

12.0 

0.75 

0 

13.0 


10 

11.0 


100 

10.5 


N.R. 

10.5 

0.50 

0 

10.5 


10 

9.5 


N.R. 

9.5 

0.25 

0 

8.0 


N.R. 

8.0 

0.10 

0 

7.0 


N.R. 

7.0 

2) r > 1.0 



(i) a = 5.0 





V 

r 

t 


2.0 

0 

116 


10 

67 


100 

15.5 


1000 

12.0 


N.R. 

12.0 

5.0 

0 

190 


10 

112 


100 

19.0 


1000 

12.0 


N.R. 

12.0 



r 

t 

Isa 

10.0 

0 

810 


10 

500 


100 

45 


1000 

12.0 


N.R. 

12.0 

50.0 

0 

20,500 


10 

12,800 


100 

1,450 


1000 

8.5* 


N.R. 

12.0 

a = 10.0 



r 

t_ 

V 

eq 

2.0 

0 

48.0 


10 

30.0 


100 

13.0 


N.R. 

12.0 

5.0 

0 

350 


10 

230 


100 

25 


1000 

12.0 


N.R. 

12.0 

10.0 

0 

1500 


10 

950 


100 

100 


1000 

12.0 


N.R. 

12.0 

50.0 

0 

37,000 


10 

25,000 


100 

2,700 


1000 

7.0* 


N.R. 

12.0 


* At t = 1000 yrs the alpha emission is greater than the beta emission 
whereas for 0, 10, and 100 yrs, it is not. 
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